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Quantum dots are small conducting devices containing up 
to several thousand electrons. We focus here on closed dots 
whose single-electron dynamics are mostly chaotic. The meso- 
scopic fluctuations of the conduction properties of such dots 
reveal the effects of one-body chaos, quantum coherence and 
electron-electron interactions.|^ 

I. QUANTUM DOTS 

Recent advances in materials science have made pos- 
sible the fabrication of quantum dots, submicron-scale 
conducting devices containing up to several thousand 
electrons 0. A 2D electron gas is created at the in- 
terface region of a semiconductor heterostructure (e.g., 
GaAs/AlGaAs) and the electrons are further confined to 
a small region by applying a voltage to metal gates, de- 
pleting the electrons under them. Insofar as the motion of 
the electrons is restricted in all three dimensions, a quan- 
tum dot may be considered a zero-dimensional system. 
The transport properties of the dot, i.e., its conductance, 
can be measured by connecting it to external leads. A 
micrograph of a quantum dot is shown in Fig. |^(a). 

At low temperatures, the electron preserves its phase 
over distances that are longer than the system's size, i.e., 
> L, where is the coherence length and L is the 
linear size of the system. Such systems are called meso- 
scopic. Elastic scatterings of the electron from impurities 
generally preserve phase coherence, while inelastic scat- 
terings, e.g., from other electrons or phonons, result in 
phase breaking 

When the mean free path £ is much smaller than L, 
transport across the dot is dominated by diffusion, and 
the system is called diffusive. In the late 1980s it became 
possible to fabricate devices with little disorder where 
£ > L. In these so-called ballistic dots, transport is dom- 
inated by scattering of the electrons from the boundaries. 
A schematic illustration of a ballistic dot is shown in Fig. 
0(b). 

In small dots (with typically less than 20 electrons), 
the confining potential is often harmonic-like, leading to 
regular dynamics of the electron and shell structure that 
can be observed in the addition spectrum (i.e., the energy 
required to add an electron to the dot). Maxima in the 
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addition spectrum are seen for numbers of electrons that 
correspond to filled (TV = 2, 6, 12) or half-filled (TV = 
4, 9, 16) valence harmonic-oscillator shells [|]. 

Dots with a large number of electrons (TV > 50 — 100) 
often have no particular symmetry, and their irregular 
shape results in single-particle dynamics that are mostly 
chaotic. For such dots the conductance and addition 
spectrum display "random" fluctuations when the shape 
of the dot or a magnetic field are varied. This is the sta- 
tistical regime, where we are interested in the statistical 
properties of the dot's spectrum and conductance when 
sampled from different shapes and magnetic fields. For 
a recent review of the statistical theory of quantum dots 
see Ref. §. 
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FIG. 1. Quantum dots: (a) a scanning electron micro- 
graph of a dot used by Folk et al. [Q. A 2D electron gas 
is formed in the interface of a GaAs/AlGaAs heterostruc- 
ture (darker area). Electrostatic potentials applied to metal- 
lic gates (lighter shade) confine the electrons to a sub-micron 
region. The shape and area of the dot can be changed by 
controlling the gate voltages Vgi and V52; (b) a schematic 
drawing of a ballistic dot attached to two leads. The elec- 
tron's trajectory scatters from the dot's boundaries several 
times before exiting. 

Many of the physical parameters of a quantum dot 
can be experimentally controlled, including its degree of 
coupling to the leads, shape, size, and number of elec- 
trons. When the dot is "open", i.e., strongly coupled to 
leads, there are generally several channels in each lead 
and the conductance fluctuates as a function of, e.g., the 
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Fermi momentum of the electron in the leads (see Fig. 
^(a)). As the point contacts are pinched off, the coupHng 
becomes weaker and a barrier is effectively formed be- 
tween the dot and the leads. In such "closed" dots, the 
charge is quantized. At low temperatures, the conduc- 
tance through a closed dot displays peaks as a function 
of gate voltage (or Fermi energy); see, e.g.. Fig. |^(b). 
Each peak represents the addition of one more electron 
into the dot. In between the peaks, the tunneling of an 
electron into the dot is blocked by the Coulomb repul- 
sion of electrons already in the dot, an effect known as 
Coulomb blockade. In this paper we discuss closed dots 
in which mesoscopic phenomena are determined by the 
interplay between single-particle chaos, quantum coher- 
ence and electron-electron interactions. 
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the confining potential written in terms of a gate voltage 
Vg and a = Cg/C, where Cg is the gate-dot capacitance. 

At low temperatures, conductance occurs by resonant 
tunneling through a single-particle level in the dot. As- 
suming energy conservation for the tunneling of the Af-th 
electron we have, Ep + £g,sXN — 1) = £g.s.(A/'), where 
Ep is the Fermi energy of the electron in the leads and 
£gs {N) is the ground state energy of a dot with M elec- 
trons. Using Eq. (|^), we find that the effective Fermi 
energy Ep = Ep 



eaVg satisfies 



Ep = E^f + [JV - 
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The conductance displays a series of peaks at values of 
Ep given by (||), with each peak corresponding to the 
tunneling of an additional electron into the dot. The 
spacings between the peaks are given by 



A2 = AEp = {Em+1 - Ej^) + eVC 



(3) 



Since the charging energy is usually much larger than the 
mean-level spacing A, the Coulomb-blockade peaks are 
almost equidistant. Coulomb blockade is illustrated in 
Fig. I 
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FIG. 2. (a) Conductance vs the electron's Fermi mo- 
mentum kp in an open dot (from Ref. Q). (b) Conduc- 
tance vs gate voltage in a closed dot displaying a series of 
Coulomb-blockade peaks (from Ref. @). 



II. TRANSPORT IN THE COULOMB 
BLOCKADE REGIME 

The simplest model for describing the Coulomb block- 
ade regime is the constant interaction (CI) model, in 
which the Coulomb energy is taken to be e^M'^j^C, 
where C is the total capacitance of the dot and M is the 
number of electrons. The Hamiltonian of the CI model 
is given by 
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Here a\\Q) is a complete set of single-particle eigenstates 
in the dot with energies E\, and J\f = "^al^ax is the 
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electron number operator in the dot. The quantity aVg is 
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FIG. 3. Schematic diagram of Coulomb blockade in closed 
dots: (a) when the gate voltage is eaVg = Ne^ jC there is a 
charging energy gap on both sides of the Fermi energy, block- 
ing the tunneling of electrons into the dot; (b) when the gate 
voltage increases to eaVg — (TV -I- l/2)e'^/C, the charging en- 
ergy for adding an electron to the dot vanishes. When the 
Fermi energy of the electron in the lead matches the first un- 
occupied single-particle state in the dot (i.e., Ep ~ Ej^+i), 
resonant tunneling of an electron into the dot occurs. From 
Ref. |. 

The Coulomb-blockade peak heights contain informa- 
tion about the wave functions. For closed dots, a typical 
level width F is small, and F ^ T even for the lowest 
electron temperatures attained in the experiments. Un- 
der this condition, the coherence between the electrons 
in the leads and the dot can be ignored, and a master- 
equation approach is feasible Q. In the CI model and 
for F ^ T <C A, the Af-th conductance peak G occurs 
through level X = J\f and is given by 
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Eq. (Q) describes a peak centered at Ep — E\ (here Ep 
is measured with respect to {N — l/2){e'^ /C)). The peak 
has a width ~ T and a height of 
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where 



3A = - 



Ta^a 
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The quantities F'^''^ describe the partial width of level A 
to decay into the left (right) lead, and F is the average 
level width. 

The partial widths can be written as Fa = |7aP, where 
7a is the partial width amplitude. In i?-matrix theory 
l^-H, 7a can be related to the wave function 5* a at the 
respective dot-lead interface 
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is the longitudinal channel momentum 
(h'^k^/2m + ?i^K^/2m = E), is the penetration fac- 
tor to tunnel through the barrier in channel c {Pc = 1 in 
the absence of barrier and Pj, *C 1 in the presence of a 
barrier), and (pc is the transverse channel wave function. 
The integral in ^ is over the dot-lead interface C. 

III. SIGNATURES OF CHAOS IN CLOSED DOTS 

In a dot where the classical dynamics of the electron 
are chaotic, the wave function fluctuations are described 
by random-matrix theory (RMT) These fluctuations 
lead to fluctuations in the conductance peaks according 
to Eqs. (|^) and (^). When time-reversal symmetry is 
conserved (i.e., there is no external magnetic field B) 
the corresponding random-matrix ensemble is the Gaus- 
sian orthogonal ensemble (GOE), while for broken time- 
reversal symmetry {B ^ 0) the appropriate ensemble is 
the Gaussian unitary ensemble (GUE). For a recent re- 
view of RMT and its applications see Ref. 0. For a 
review of the random-matrix theory of quantum trans- 
port (including applications to open dots) see Ref. Jill . 

To quantify the conductance peak statistics using 
RMT, we express the partial width amplitude (^ as 
a scalar product of the resonance wavefunction tpx = 
('(/'Ai, '0A2, • ■ ■) and the channel wavefunction (f)c = 
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where we expanded the wavefunction 'I'a = ^ ^^jPj 

j 

a fixed basis pj in the dot and defined the channel vector 



scalar product in (|^) is defined over the dot-lead inter- 
face and is different from the usual scalar product in the 
Hilbert space of the dot. 



A. Peak heights distributions 
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FIG. 4. Conductance peak-height distributions 
in Coulomb-blockade quantum dots: (a,b) measured distri- 
butions by Chang et al. (l3| for (a) B = and (b) B / 0. 
The effective size of the dot is 0.25 x 0.25 /im and T — 75 
mK (corresponding to TV ~ 100 electrons and T/A ~ 0.15). 
The solid lines are the RMT predictions of Jalabert, Stone 
and Alhassid (c,d) measured distributions of Folk et al. 

compared with the RMT distributions (solid lines). The 
dots used have areas of 0.32 and 0.47 /im^ and T = 70 ± 20 
mK {Af - 1000 and T/A ~ 0.3 - 0.5). Adapted from Refs. 
H and H 
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The conductance peak distributions for one-channel 
symmetric leads (A = 1; r' = f '') were derived by Jal- 
abert, Stone and Alhassid using RMT ||]. Similar results 
were obtained in the supersymmetry method [ p^ . Using 
Eqs. (|) and (|), and assuming ^\ to be an RMT eigen- 
vector, we can calculate the distributions of the dimcn- 
sionless conductance peak heights g. These distributions 
are universal and depend only on the space-time symme- 
tries iO: 
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where Kq and Ki are modified Bessel functions. The 
distributions were measured independently by Chang et 
al. jl^ and by Folk et al. and found to agree with 
the theoretical predictions for both conserved and broken 
time-reversal symmetries. Fig. || shows a comparison 
between the theoretical and experimental distributions 
in both experiments. The peak height distributions in 
closed dots with several possibly correlated and inequiv- 
alent channels in each lead were derived in Refs. Ifl5,[l6i. 



B. Weak localization 

Another signature of chaos is the weak localization 
effect. This well-known quantum interference effect 
was already observed in macroscopic disorder conduc- 
tors jr^ . When time-reversal symmetry is conserved, 
time-reversed orbits contribute coherently to enhance the 
return probability. The average conductance is then 
smaller in the absence of magnetic field. 




FIG. 5. Weak localization in Coulomb-blockade dots. The 
analytical prediction (^o|) for the average conductance g/g'^^^ 
vs a time-reversal symmetry breaking parameter ^ (solid line) 
[ po[ is compared with the experimental results of Ref. ^] 
(solid circles). In the experiment — B/B^r is a scaled mag- 
netic field {Bcr ~ 6 mT), and g'^^^ is measured for B away 
from zero. 



In Coulomb-blockade quantum dots, the average con- 
ductance peak height can be calculated as a function of a 
time-reversal symmetry breaking parameter ^. In RMT, 



the statistics of the crossover regime between GOE and 
CUE can be described by the Mehta-Pandey ensemble 



H = S + iaA 
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where S and A are, respectively, symmetric and an- 
tisymmetric real matrices and a is a real parameter. 
The matrices S and A are uncorrelated and chosen 
from Gaussian ensembles of the same variance. The 
transition parameter C is given by a typical symmetry- 
breaking matrix element measured in units of A [ p^ : 
( = (inbreak) A = ay/N/n (where N is the dimen- 
sion of the random matrix). If the time-reversal symme- 
try is broken by a magnetic field B, then C = B/Bcr, 
where Bcr is the crossover field. 

Using the ensemble (^) to describe the statistics of the 
partial amplitude (^, we can find the average conduc- 
tance peak height in a dot with single-channel symmetric 
leads iol 



m = \ 
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Eigenvectors in the crossover ensemble have complex 
components whose distribution in the complex plane is 
characterized by a "shape" parameter t. The function 
P({t) in ( pO| ) describes the probability of an eigenvector 
to have a certain shape t and is given by [pTl-^ 

Pc{t)=7T 




where (fii {() = J e 



The average dimensionless conductance 5(C) has a dip 
at C = (see solid line in Fig. |^), describing a weak 
localization effect. The theoretical results compare well 
with recent experimental results vs magnetic field p^ ] 
(solid circles in Fig. |5|) once the magnetic field is scaled 
by Bcr ~ 6 mT. 

The crossover field can be estimated semiclassically. 
Time-reversal symmetry is fully broken for field Bcr 
where the rms of the phase difference between an or- 
bit and its time-reversed partner is ~ 27r. This phase is 
proportional to the area enclosed by the electron's trajec- 
tory. In a chaotic dot, area accumulation is diffusive, and 
the accumulated area's rms behaves as the squared-root 
of the elapsed time . In an open dot the relevant time 
is the escape time Tcscapo, but in a closed dot this time 
is replaced by the Heisenberg time th = h/ A. A more 
quantitative derivation (for a dot with area A) gives 
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(12) 

where is the ergodic time (roughly the time to cross the 
dot), and — Et / IS. (with Et = fi/rc) is the baUistic 
Thouless conductance. The factor k is a non-universal 
geometrical factor. 

C. Parametric correlations 

Another signature of quantum chaos are the meso- 
scopic fluctuations of a given conductance peak height 
as a function of an external parameter, e.g., the shape 
of the dot or a magnetic field. These fluctuations can be 
described in the framework of Gaussian processes (GP) 
p6t which generalize the random-matrix ensembles to 
random-matrix processes. A Gaussian process H{x) of 
a given symmetry class (3 is characterized by its first two 
moments 



H,^{x)Hki{x') = ^f{x~x')g^l, , 



where the coefficients g\^\i are defined by 
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The Gaussian process (13) constitutes a Gaussian ensem- 
ble for each value of the parameter x. 

A GP is characterized by the short distance behavior 
of its correlation function f{x — x') ^ \ — k\x — x'Y' . 
A differentiable GP is obtained for rj — 2 |2^, and cor- 
responds to the usual situation where the Hamiltonian 
depends analytically on the parameter. Simons and Alt- 
shuler showed that parametric correlations in dis- 
ordered or chaotic systems become universal once the 
parameter x is scaled by the rms of the level velocity pq] 
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Here is the ith energy level measured in units of the 
mean-level spacing A. The conductance peak correla- 
tor is 

defined by Cg{x - x') = 6G{x)6G{x')/[a{G{x))c 7{Gix'))l 
where 6Gix) = G(x) - G{x) and <y^{G{x)) = WH^W 
(here G{x) is the conductance peak height at a value x 
of the parameter). The universal correlator Cg was cal- 
culated in Ref. using the simple GP @|6| H{x) = 
Hi cosx + H2 since, where H2 are uncorrelated ma- 
trices that belong to the appropriate Gaussian ensemble. 
The GUE correlator can be well fitted to the square of a 
Lorentzian. If the parameter is a magnetic field B then 
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where Be is the correlation field. Fig. |6| shows the mea- 
sured conductance peak correlator (solid diamonds) in 
comparison with the theoretical prediction ( |T^ ) (solid 
line), where the parameter Be is fitted. RMT describes 
the correct shape of the correlator for B < 10 mT. A 
single-particle semiclassical estimate of the correlation 
field Be gives an expression similar to (^2|). For a sta- 
dium in a uniform magnetic field, the correlation flux 
is ~ 0.3 $0 23|, which is below the experimen- 
tal value of « 0.8 $o- This indicates that the single- 
particle picture is inadequate for estimating the correla- 
tion field. Numerical simulations in small disordered dots 
with Coulomb interactions find that the correlation field 
increases with the interaction strength [^^ . 
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FIG. 6. The conductance peak height correlator Cg vs 
magnetic field AB in closed dots. The experimental results of 
Ref. (solid diamonds) are compared with the theoretical 
prediction (^ of Ref. Ill (solid line). The inset shows the 
measured peak height of a single Coulomb-blockade peak vs 
magnetic field B. Adapted from Ref. [IfJ. 



D. Peak-to-peak correlations 

The conductance peak at finite temperature (T 
can be calculated in the master-equations approach 



A) 



G{T,Ef) 



where g = Y.xWx{T,EF)gx (17) 



is the dimensionless conductance expressed as a ther- 
mal average over the level conductances gx = 
2f -^T[Tl/{T{ + Tl). For T, A < e^/C, the thermal 
weights for the M-th conductance peak are given by 



Ef 



(18) 



where P_\f is the probability that the dot has J\f electrons, 
and {n\) is the canonical occupation of a level A. 

The finite-temperature peak height statistics were cal- 
culated in Ref. |Q assuming the level conductance and 
energy levels satisfy RMT. Of particular interest is the 
peak-to-peak correlator 



5 



c{n) = 5GM5GM+n/{5G^rY , 



(19) 



where 5Gj^ — Gj^ — Gj^f is the fluctuation of the M- 
th conductance peak around its average. The number of 
correlated peaks Uc is the FWHM of ( [l^ . The calculated 
temperature dependence of iic is shown in Fig. ^c) (solid 
line). 

The experimental results shown in Fig. ^ demon- 
strate that ric saturates with temperature, especially for 
the small dots. This behavior is contrary to the linear 
dependence predicted by RMT. An explanation of this 
saturation effect is discussed in Sec. [VA. 
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FIG. 7. Finite-temperature peak-to-peak correlations (ex- 
periment). The peak-to-peak correlator c(n) is shown at dif- 
ferent temperatures for (a) a larger dot with A = 20 /leV and 
(b) a smaller dot with A = 38 /leV. (c) The number of cor- 
related peaks ric (defined as the FWHM of c(n)) vs T/A for 
three dots. The solid line is the RMT result for an unchanged 
spectrum. From Ref. M . 



IV. INTERACTION EFFECTS 

While the CI-plus-RMT model can explain some of the 
observed statistical properties of closed dots, several ex- 
perimental observations suggest that electron-electron in- 
teractions beyond the charging energy are important: (i) 
The measured peak-spacing distribution [^ -37| does not 
have the Wigner-Dyson shape expected in the CI model 
(see Eq. (|^)). This signature will be discussed in detail 



in Sec. IV B. (ii) The measured correlation flux is larger 
than the semiclassical estimate (see Sec. Ill C] ). (iii) Cor- 
relations between the addition and excitation spectra are 
seen only for a small number of added electrons . (iv) 
The peak-to-peak correlations saturate with increasing 
temperature, contrary to the results of the CI-plus-RMT 
model. This effect is explained in Section [VA. 



A. Spectral scrambling 

The best way to include interaction effects while re- 
taining a single-particle framework is the mean-field ap- 
proach, e.g., the Hartree-Fock approximation. In this 
approach, the self-consistent single-particle energy levels 
change ("scramble") when an electron is added to the 
dot. This scrambling can be modeled in terms of a dis- 
crete parametric random matrix H(x/^), where de- 
scribes the parameter of the dot with J\f electrons. We 
describe H{xj\f) as a discrete Gaussian process, and em- 
bed it in a continuous process H{x). We assume that 
the scaled parametric change 6x = xa^+i — xj\f upon the 
addition of one electron is independent of Af. 

The parameter Sx determines the degree of spectral 
scrambling. Fig. ^(a)(b) shows the correlators c(n) at 
various temperatures for both an unchanged spectrum 
{6x — 0) and for a spectrum that scrambles {6x — 0.5). 
Fig. 8(c) shows that when 6x is larger (i.e., the spec- 
trum scrambles faster), the number of correlated peaks 
saturates at smaller values of T/A 

Spectral scrambling is an interaction effect [Q. A 
microscopic estimate of Sx can be done in Koopmans' 
limit pO[ |, i.e., assuming the single-particle wave func- 
tions do not change with the addition of electrons (and 
only the spectrum scrambles). In this limit the change 
SEi of an energy level upon the addition of one electron is 
given by a diagonal matrix element SEi ~ '^i.A/'+i: where 
Vaj = Vajia-y- Usiug an RPA screened interaction |4l| ] or 
a short-range dressed interaction v ^ XAAS{r — r'), 
we have 



(20) 



where gr is the Thouless conductance. In the finite dot, 
there could be an additional contribution due to excess 
negative charge on the boundaries M], leading to 



aHSE.)r.^^ 
P 9T 



(21) 



In the parametric approach a^(SEi) w (Sx)'^, and by 
comparing with (|2^) or we can find Sx. The number 
of added electrons m required for complete scrambling 
of the spectrum is then determined from nicSx ~ 1. We 
obtain MM 
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where (^2b|) holds in the presence of surface charge 
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FIG. 8. Finite-temperature peak-to-peak correlations and 
spectral scrambling (theory). The calculated correlators c(n) 
are shown at several temperatures for (a) a dot with an un- 
changed spectrum {5x = 0) , and (b) a dot where the spectrum 
scrambles [Sx — 0.5). (c) The number of correlated peaks 
ric vs T/A for dots with different scrambling parameters 5x. 
Notice the earlier saturation of as the spectrum scrambles 
faster (i.e., when 5x is larger). From Ref. [M. 



B. Peak spacing statistics 

One of the main signatures of interactions is seen in the 
peak spacing distribution. The spacing between peaks in 
an interacting dot is given by 

A2 = f g.s. £g.s. (AA - 1) - 2£g,s. (M) , (23) 

where £g.s. (A/") is the ground state energy of the dot 
with A/" electrons. In the CI model, A2 reduces to (|^), 
and a shifted Wigner-Dyson distribution is expected for 
A2. Experimentally, the distribution is Gaussian-like (for 
semiconductor dots with ~ 1 — 2) and has a larger 
width than the Wigner-Dyson distribution ||35|-|37||. An 
example of the measured peak spacing statistics is shown 
in Fig. |. 

An estimate of the fluctuations of A2 can be obtained 
in Koopmans' limit p|] where the addition energy is 

given by E^J^'^ (-^l"^'' is the i-th single-particle level of 
a dot with Af electrons) . The peak spacing can then be 
written as 
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(24) 



The first term on the r.h.s. of (gj) is the usual level spac- 
ing in a dot with a fixed number of electrons {M + 1) and 
has an rms of order A. The second term A.E^ represents 
the change of a given energy level when an electron is 
added and its standard deviation can be estimated from 
® or dll). 
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FIG. 9. Measured peak-spacing statistics in quan- 
tum dots: (a) a sequence of peak-spacing fluctuations 
Us = {AVg ~ 'AV^)/AVg vs. Vg iov B = 30 mT. The dashed 
lines show the standard deviation of RMT. (b) The mea- 
sured peak-spacing distribution P{vg) (shaded histograms) 
for B = and T/A ~ 0.8 is compared with a Gaussian fit 
(solid line). From Ref. [||. 

Gaussian-like peak spacing distributions were ex- 
plained in exact numerical diagonalization of a small An- 
derson model with Coulomb interactions [ ^5p^ . Their 
width increases with the gas constant Tg. Scrambling of 
the single-particle spectrum (see Section IVA) can also 
lead to Gaussian distributions if 5x is sufficiently large 

El- 

Nearly-Gaussian dis- 
tributions were observed in Hartree-Fock calculations at 
larger values of |^6|-Q . The distribution of a diagonal 
interaction matrix element is found to be approximately 
Gaussian. Using (|2^ ) (for small values of Ts), we can de- 
scribe the peak spacing distribution as a convolution of 
the Wigner-Dyson distribution with a Gaussian distribu- 
tion. 

The numerical investigations of small disordered dots 
with interactions demonstrate the need to go beyond the 
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simple CI model. An interesting question is whether it is 
possible to describe the RMT-like behavior of the peak 
height statistics and the Gaussian-like distribution of the 
peak spacing distribution within a single random-matrix 
model. 

C. Random interaction matrix model 

RMT is not restricted to single-particle systems. It was 
successfully applied to strongly interacting systems, e.g., 
the compound nucleus at finite excitations. However, in 
the linear conductance regime of quantum dots we are in- 
terested in the statistical properties of the ground state 
of the system as the number of electrons changes. Since 
RMT does not make explicit reference to interactions or 
to particle number, it is necessary to use a model that 
contains interactions explicitly. A two-body random- 
interaction model was introduced in nuclear physics in 
the early 1970s [ ^9|j5^ . It was used, together with a 
random single-particle spectrum, to study thermalization 
pl| and the onset of chaos in interacting many-body sys- 
tems . Since its one-body part has Poissonian statis- 
tics it is not suitable for studying chaotic dots. A random 
interaction matrix model (RIMM) for dots with chaotic 
single-particle dynamics was recently introduced to study 
the interplay between one-body chaos and interactions 

The RIMM describes an ensemble of interacting Hamil- 
tonians in a fixed basis of m single-particle states \i) — 

«llo) 

^ ^Yl ^^i^-l^j + ^tjkA^^-iak ■ (25) 

ij ijkl 

The one-body elements hij are chosen from the appropri- 
ate Gaussian random-matrix ensemble, i.e., GOE (GUE) 
for conserved (broken) time-reversal symmetry of the 
one-body dynamics. The anti-symmetrized two-body 
matrix elements uf^.^-i = Uij-^ki —Uij-^ik form a GOE in the 
two-particle space (the two-body interaction is assumed 
to conserve time-reversal symmetry irrespective of the 
symmetry of the one-body Hamiltonian) . The variance 
of the diagonal (off-diagonal) interaction matrix elements 
is {U'^/2). The two-body interaction can include a 
non-vanishing average part u that is invariant under or- 
thogonal transformations of the single-particle basis. The 
only such invariant for spinless electrons is the charging 
energy e^Af^ /2C, which is a constant and thus does not 
affect the statistical fluctuations. If the spin degrees of 
freedom are included in the RIMM then another possible 
invariant is the exchange interaction |54j. It is 

important to note that in a physical model of the dot, the 
two-body interaction in a given basis is fixed. The intro- 
duction of a fluctuating two-body interaction is done in 



the spirit of RMT to describe generic effects that do not 
depend on the specific interaction. 

The RIMM was used to calculate the peak spacing dis- 
tribution from Eq. (Eq), where the ground state energies 
of the Hamiltonian (Ea) are calculated for A/" — 1, A/" and 
N +1 electrons. The peak spacing distribution describes 
a crossover from a Wigner-Dyson distribution at J7 = to 
a Gaussian distribution as J7/A increases; see Fig. px|(a) 
p3[ . The width of the distributions increases vs C//A 
(see Fig. |l^(b)). The distributions are well described 
by a convolution of a Wigner-Dyson distribution with a 
Gaussian distribution. For small J7/A, the width of this 
Gaussian is just the standard deviation of a diagonal in- 
teraction matrix element. 

The inset of Fig. |l^(b) shows the calculated ratio 
crGOE(A2)/o'GUE(A2) VS C//A. The experimental values 
~ 1.1 — 1.3 are consistent with theory. 



0.8 




12 3 4 

U/A 



FIG. 10. Peak spacing statistics in the RIMM: (a) B = 
peak-spacing distribution for several values of C//A: 0.35 
(open squares), 0.7 (solid diamonds), 1.1 (open triangles) and 
1.8 (solid triangles). The dashed line is the Wigner-Dyson dis- 
tribution and the solid line is a Gaussian fit to the U /A = 1.8 
distribution, (b) The standard deviation a(A2) vs C//A for 
GOE (solid line) and GUE (dashed line) one-body statistics. 
The inset is the ratio (tgoe(A2)/(tgue(A2) vs f//A. All the 
results are shown for m = 12 and A" = 4. Adapted from Ref. 

il. 

To calculate the conductance peak height we use Eq. 
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but now the partial width of the ground state of the 
AA-electron dot to decay into the ground state of the dot 
with M — 1 electrons plus an electron in the respective 
lead is given by: 



(26) 



where ip\r) is the creation operator of an electron at the 
point contact r, and ^g,s.{-^) is the ground state wave- 
function of the A/'-electron dot. For a GOE one-body 
statistics we find that the partial-width distribution is 
a Porter-Thomas distribution irrespective of U /A (see 
Fig. 11(a)). A similar insensitivity to the interaction 



strength was found in numerical simulations of an An- 
derson model with Coulomb interaction s [^ . However, 
for a GUE single-particle statistics (Fig. |ll|(b)) we find a 
crossover from the GUE Porter-Thomas distribution (at 
U=0) to the GOE Porter-Thomas distribution (at large 
U/A). Equivalently, the GOE GUE transition due to 
an external magnetic field is incomplete because of the 
competing GOE symmetry of the two-body interaction. 
The crossover width distributions are well described by 
those of the Mehta-Pandey ensemble (||) with a transition 
parameter that is a monotonically decreasing function 
of C//A II . 



FIG. 11. Width and conductance peak height statistics 
in the RIMM: (a) Width distributions P(lnf) vs. Inf are 
shown for GOE one-body statistics and U/A = (circles), 
2.4 (squares), and 4 (diamonds). The solid and dashed lines 
are the GOE and GUE Porter- Thomas distributions, respec- 
tively, (b) Same as in (a) but for GUE one-body statistics. 
The short-dashed lines are analytic width distributions de- 
rived from the crossover ensemble (^. All the results are 
shown for m = 12 and J\f = 4. The insets in (a) and (b) are 
the corresponding peak-height distributions P{g) in a log-log 
scale. Adapted from Ref. |p5|. 

We note that the curves of the width cr(A2) and C vs 
U/A depend on both m (number of single-particle states) 
and J\f (number of electrons) , but they become universal 
once U/A is scaled by a constant (that depends on m and 
Af). For U/A ^ 0.7—1.5, the peak spacing distribution is 
already Gaussian-like, while the peak height statistics (in 
the presence of a time-reversal symmetry breaking field) 
is still close to the GUE prediction. In the RIMM, U/A 
is a free parameter and physical values can be determined 
by comparing its results to specific models. For a small 
(4 X 5) Anderson model with Coulomb interactions of 
strength Uc, we find that the range C//A ~ 0.7 — 1.5 
corresponds to Uc ^ 2 — 5 l53,p3| . 




V. CONCLUSIONS 

The mesoscopic fluctuations in closed dots are affected 
by both one-body chaos and electron-electron interac- 
tions. Signatures of chaos include the RMT-like conduc- 
tance peak height distributions, the weak localization ef- 
fect, and the line shape of the parametric peak- height 
correlator. Interaction effects include the Gaussian-like 
shape of the peak-spacing distribution, a correlation 
magnetic field that is larger than its single-particle es- 
timate, and the saturation of peak-to-peak correlations 
with temperature. Some of these interaction effects can 
be described by a random interaction matrix model. 
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